Gravitational Radiation from Preheating with Many Fields 
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Parametric resonances provide a mechanism by which particles can be created just after inflation. 
Thus far, attention has focused on a single or many inflaton fields coupled to a single scalar field. 
However, generically we expect the inflaton to couple to many other relativistic degrees of freedom 
present in the early universe. Using simulations in an expanding Friedmann-Lemaitre-Robertson- 
Walker spacetime, in this paper we show how preheating is affected by the addition of multiple fields 
coupled to the inflaton. We focus our attention on gravitational wave production-an important po- 
tential observational signature of the preheating stage. We find that preheating and its gravitational 
wave signature is robust to the coupling of the inflaton to more matter fields. 



I. INTRODUCTION 

After inflation the universe is cold and far from thermal 
equilibrium. Early ideas for how to reheat the universe 
after inflation involved standard model particle produc- 
tion via oscillations of the inflaton about the minimum 
of its potential The coupling between the inflaton 

and matter fields, however, can lead to parametric ampli- 
fication of the quantum fluctuations of the matter fields. 
Thus, the reheating process can be preceded by a stage of 
exponential particle production which has become known 
as preheating [5H7|. 

Inflation produces a stochastic background of gravi- 
tational waves through the amplification of primordial 
quantum fluctuations @, Q . It turns out that preheating 
may also lead to the production of a gravitational wave 
background. This realization was first made by Khleb- 
nikov and Tkachev (To| in a quartic inflation model. A 
similar estimate by Garcia-Bellido somewhat later [Tl| . 
showed that in hybrid inflation the frequency of the peak 
in the spectrum could be brought down to about 1 kHz, 
the frequency range where where current earth-based in- 
tcrfcromctric gravitational-wave detectors operate. 

This problem has received much attention recently [l2| - 
[20j . Most of this work uses LatticeEasy [2l| to nu- 
merically evolve the inflaton and other scalar fields in 
Friedmann-Lemaitre-Robertson- Walker (FLRW) space- 
times [21], but different schemes to compute the metric 
perturbation and gravitational wave spectrum. There is 
now general agreement on the gravitational wave signal 
predicted by the preheating mechanism [l9[. All of this 
work used toy models with one inflaton field and one 
(scalar) matter field. More recently there has been some 
progress in our understanding of preheating with many 
inflaton fields [22Tj24l |. a possibility suggested by string- 
theory-based inflation scenarios. 

Generally, we expect the inflaton to couple to many of 
the degrees of freedom present in the early universe. Up 
until now, the effect of these extra fields on the resonance 
and the gravitational wave signature of preheating has 
not been explored. Using simulations in an expanding 



FLRW spacetime, in this paper we show how preheat- 
ing is affected by the addition of multiple fields coupled 
to the inflaton, focusing our attention on gravitational 
wave production. We briefly discuss the possibility of 
determining the number of fields the inflaton is coupled 
to from the gravitational wave signature. 

In Sec. [II] we review the dynamics of preheating and 
outline our computational strategy to calculate the grav- 
itational wave spectrum. In Sec. lIIII we discuss our results 
for quartic and quadratic models of inflation, for the case 
when the inflaton is coupled to between 1 and 32 matter 
fields. Throughout we compare our results to the well- 
known case where the inflaton is coupled to a single field. 
We conclude in Sec. HVl 



II. MODEL AND COMPUTATIONAL METHOD 

We will consider a simple extension of the usual pre- 
heating model [5[ in which the inflaton, <p, is coupled to 
a set of Af scalar fields, Xa, according to the Lagrangian 



(i) 

where V((f),Xa) contains both the inflationary potential, 
V((f>) and the interactions between the inflaton and the 
scalars, 



.V 



(2) 



For simplicity we will assume that the coupling between 
the matter fields and the inflaton, g 2 , is the same for all 
fields. The fields evolve in a FLRW spacetime, which 
takes the form 

ds 2 = dt 2 - a 2 (t) [dx 2 + dy 2 + dz 2 ] (3) 
= a 2 (r) [dr 2 - dx 2 - dy 2 - dz 2 ] , (4) 
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in cosmological and conformal coordinates respectively. 
In this paper we will consider both quadratic, V((j)) = 
m 2 cj> 2 /2, and quartic, V(<j>) — Xcj) 4 /^, inflationary models. 
The fields obey the usual Klein-Gordon equations, 



3H(f> 



v 2 ^ dv(<i>) 

a 2 d<fi 



AT 



J2g 2 txi = o, (5) 



and 



X a + 3HXa ~ 5" + g 2 2 Xa = 0. 



(6) 



To show that we can achieve parametric resonance in this 
model, it's easiest to recast Eq. (j6|) in terms of the mode 
equations of Xa, 

i + 3H Xa + + ffV ) Xa = 0, (7) 



where we use the Fourier convention, 

/(M)= / d 3 xf(x,t)e 2 ^ s . 



(8) 



Before we extend our analysis to many fields, we will 
briefly review the phenomenon of parametric resonance. 
In both quadratic and quartic inflation, the end of in- 
flation is marked by coherent oscillations of the 4> field. 
At the end of quadratic inflation, temporarily ignoring 
the expansion of the universe, the inflation will oscillate 
sinusoidally, 



4>(x,t) — <j>sin(mt). 



(9) 



In the case of a single coupled scalar, y, the mode equa- 
tions Eq. (O can be cast into a Mathieu equation j25[, 



x" + (A fe -2gcos(2z))x = 



using the coordinate transformation, 



q = 



g 2 * 2 

4m 2 



g 2 ® 2 
2m 2 



mt, 



(10) 



(11) 



where primes denote differentiation with respect to z. 
Solutions to Eq. (fTOj) can be characterized as either stable 
(oscillatory) or unstable (exponential) depending on the 
choice of parameters 1 , Af. and q. Preheating, in the form 
of parametric resonance, occurs when particular modes of 
X are excited, producing a universe whose content is out 
of thermal equilibrium 0, 0] ■ I n practice, the inclusion of 
the expansion of the universe complicates the coordinate 
transformation Eq. (jll[) by introducingtime-dependent 
Mathieu parameters, Ak(t) and q(t) 



In the case of quartic inflation, the time dependence of 
<p is not purely sinusoidal; rather, the field solution can 
be written in terms of an elliptic cosine 



cx cn \ x 



1 

71 



(12) 



where x cx \Jt [26[. Substituting this into the single-field 
version of Eq. ([7]) yields a Lame equation. The dynamics 
of this can be understood as an approximate Mathieu 
equation since the elliptical cosine can be decomposed 
into sinusoidal parts, each of which have an associated 
Mathieu equation |26j |. 

This process is non-thermal and generates large 
anisotropics in the stress-energy tensor. Khlebnikov 
and Tkachev predicted that this anisotropic stress-energy 
should be an efficient generator of stochastic gravita- 
tional waves fl0| . The problem remained largely ignored 
for several years until Easther and Lim confirmed the ef- 
fect [IH, in work that sparked a resurgence of interest 
and activity in the field |15l420| |. 

In principle, our understanding of parametric reso- 
nance can be trivially extended to the case of many scalar 
fields, Xa- Since we are taking a single degree of freedom 
for the inflaton, the behavior of <p(t) in our model Eq. ([1} , 
at least initially, should be identical to the single-scalar 
case. The only difference being that there will be Af iden- 
tical Mathieu/Lame equations, parametrically sourcing 
inhomogeneities in each of the Xa- 

We do not expect trivial behavior in the dynamics of 
multi-scalar preheating. In the case of a single scalar, 
parametric resonance continues until the backreaction 
becomes important and the inhomogeneities in x cause 
the inflaton to loose coherence via the interaction poten- 
tial V{<j>, x). In the case of many scalars, we expect this 
effect to be enhanced. 

Although the process is highly nonlinear, qualitatively 
we expect the introduction of many scalar fields Xa 

to 

have an effect on the efficiency of preheating either by 
shortening the process or limiting the amplification of the 
modes of any particular Xa- The open questions are: (1) 
how different are the nonlinear dynamics of preheating 
in the case of many-scalars, and (2) what effect will this 
have on the generation of gravitational radiation. 

To address these questions, we have performed a se- 
ries of numerical simulations with Af = 1, 2, 4, 8, 16, and 
32 matter fields Xa- We use the publicly available early 
universe field evolution package LatticeEasy to per- 
form our simulations. We couple LatticeEasy to a code 
that evolves the metric perturbation using the methods 

of [HG1. 

Tensor perturbations to the metric, hij, can be written 
in the synchronous gauge 



ds 2 



dt 2 - a 2 (t) [6ij + hij] dx l dxK 



(13) 



For this choice of gauge perturbations are purely spatial 
and obey the transverse-traceless conditions 



1 See, for example, [TJ§ for stability diagrams. 



and hi 



0. 



(14) 
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FIG. 1: The field variances, <f> (blue) and \ (red), as a function 
of conformal time, r, for the M — 1 case. In this simulation 
the inflaton potential is quartic with A = 10 -14 and <? 2 /A = 
120. 



The metric perturbations obey sourced Klein-Gordon 
equations 



hij + 3Hhi 



16ttS 



TT 



(15) 



where the source term, Sj 3 T is the transverse-traceless 
part of the anisotropic stress tensor, 



Si 



Vij 



T. 



(16) 



We work in natural units where the speed of light and 
Newton's gravitational constant are set to unity, c = G = 
1. 

After we specify our model Eq. (JIJ, LatticeEasy 
evolves the fields and the scale factor. From this, we 
calculate the momentum-space anisotropic stress tensor, 
Eq. (|16p. which we use to source the evolution of the 
metric perturbation h^. 

We compute the present-day power spectrum of 
stochastic gravitational radiation coming from preheat- 
ing as follows. The stress energy tensor associated with 
the metric perturbations is [27], 



1 



(17) 



where the brackets denote a spatial average. The energy 
density in gravitational waves is given by 



Pgw 



— (hij nh%) = —Y (hi) , 
32-7T X ' 1 V>m-t-~l\ 



32tt 



(18) 



where V is the comoving volume over which the spatial 
average is being performed. We can then write 



dp. 



dlnfc 32 



(20) 



which is related to the present-day power spectrum 
via [Hill, 



( \ 1/3 1 d 

\9e I Ptot.e 



din k 



(21) 



where the and e subscripts denote quantities today and 
the end of our simulations respectively, h is a dimension- 
less factor that absorbs the uncertainty in the present 
value of the Hubble parameter, fi ra d,o is the current frac- 
tion of the energy density in the form of radiation, and 
Ptot.o is the total energy density at the end of our simula- 
tions. The ratio, go/g e , is the number of degrees of free- 
dom today to the number of degrees of freedom at mat- 
ter/radiation equality. We approximate ga/g c = 1/100. 

Our simulations run on a 3-dimensional lattice with 
128 points along each direction. 



III. RESULTS 

In both massive and massless models, the period im- 
mediately following inflation is characterized by the co- 
herent oscillation of the inflaton, (f>, about its minimum. 
Since this is the dominant mode, during this period the 
inflaton acts as a homogeneous field $(t). As such, the 
Xa fields each obey identical equations of motion, 



Xa + 3Hxa - 



+ g^ 2 (t) Xa =0. 



(22) 



We expect the Xa fields to undergo resonance during this 
stage, each field acting identically. 

Our expectations for the dependence on Af for the 
gravitational wave spectrum during this period can be 
sketched out with the following simple argument. The 
metric perturbation is given by 

M*. k )« f dt'G(t,t')T^ T (t',k), (23) 



where G(t,t') is the Green's function and t is the appro- 
priate time coordinate, either conformal or co-moving, 
and the overdot represents a time derivative. The gravi- 
tational wave spectrum is 



where = (1,0,0,0) in the background metric defined 
in Eq. (ITB")) and the overdot denotes a time derivative. 
Using Parseval's theorem (see [ll|) we can compute the 
power spectrum at the end of the simulation, 



1 1 

32^U 



E 



hij(t, k) 



(19) 



i,3 



gw ~ ^ dlnfc 



oc YjJ dQ,dt'dt"G(t,t')G{t,t") 

i,3 



x^T(i',k)7- TT (;",k). (24) 
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For each direction SI in the integral above, we can perform 
a rotation to a frame where the gravitational wave is 
traveling along the ^-direction. In this case the square of 
the stress-energy tensor has a particularly simple form 

J2 nr'7-' IT = (\(T XX ~ T yy )/2\f + (\T xy \f . (25) 

Every component of the stress energy tensor is a sum of 
contributions from each field, 

T H = ^2 ^iXaVjXa ~ 9ij [VfiXaV^Xa + V(<P, Xa)} ■ 
a=l 

(26) 

Since the field values are uncorrelated we expect the sum 
to grow like VAf, and hence Eq. (|2l)|) and the gravita- 
tional wave spectrum should scale like Af. 

After this initial period, the evolution is highly nonlin- 
ear and we resort to numerical simulations to study it. 
In the following two subsections we investigate the cases 
of quartic and quadratic inflation. 



A. Quartic inflation 

We start by examining massless inflation. The infla- 
tionary potential is 

V(j) = (27) 

where we take the value of the self-coupling, A = 1CP 14 , 
and g/X = 120. These par ticular values are chosen to 
make comparisons with [l6l Il7l . [l9l simple and also to 
satisfy the COBE normalization [28j . 

The dynamics observed in our simulations can be de- 
scribed in three stages. During the first stage the inflaton 
is oscillating coherently at the bottom of its potential, 
inducing parametric resonance in each matter field inde- 
pendently. The gravitational wave spectrum at this point 
in the simulations is shown in Fig. [3] At both high and 
low frequencies the expected scaling of the gravitational 
wave spectrum proportional to Af is observed. Overall 
the gravitational wave production at this stage is not 
substantial. The peaks of the spectra will rise roughly 
15 orders of magnitude before the end of the preheating 
process. 

The next stage of the evolution is highly nonlinear. It 
is marked by the amplification of perturbations of the 
inflaton driven by the matter fields. The variances of the 
inflaton fields from each of the simulations are shown in 
Fig. El The simulations with more matter fields create 
more backreaction on the inflaton, causing the largest 
fluctuations in the inflaton in the Af — 32 simulation, 
followed by the Af = 16 simulation, and so on. Gravi- 
tational wave production at this stage is significant, as 
finally summarized in Fig. Ul but the peaks of the spec- 
tra are still approximately five orders of magnitude away 
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FIG. 2: The present day gravitational wave spectrum pro- 
duced during the stage when the <j> field is coherently oscil- 
lating. The color coding is: Af — 1 (red), Af = 2 (orange), 
Af — 4 (green), Af = 8 (blue), Af = 16 (purple), Af = 32 
(black). In these simulations the inflaton potential is quartic 
with A = 10" 14 and g 2 /\ = 120 
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FIG. 3: The variances of <j> for each simulation as a function 
of conformal time, r, during the stage of coherent oscillations 
of <j>. The color coding is: Af — 1 (red), Af = 2 (orange), 
Af = 4 (green), Af = 8 (blue), Af = 16 (purple), Af = 32 
(black). In these simulations the inflaton potential is quartic 
with A = 10" 14 and g 2 /\ = 120. 



from their final values. At the same same time, the in- 
creased backreation decreases the efficiency of the reso- 
nance and the result is that the inflaton resonance ends 
sooner in the simulations with more matter fields. The 
effect on the matter fields is shown in Fig. [5J where the ef- 
fect of the different Af becomes apparent in the evolution 
of the matter fields. 

The nonlinearity of this stage makes it difficult to pro- 
duce an analytic argument to understand the details of 
the A/"-dependence. In the final stage, the fields settle 
and, in this model of inflation, thermalize. 

The final gravitational wave spectra are shown in 
Fig. |6l Note that during this stage the higher frequency 
modes are excited more in the simulations with fewer 
fields and this leads to a crossing over over of the spec- 
tra. 
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FIG. 4: The present day gravitational wave spectrum pro- 
duced during the stage when the inhomogeneities of 4> axe 
growing, and the \ fields are still resonating. The color cod- 
ing is: TV = 1 (red), Af — 2 (orange), Af = 4 (green), Af = 8 
(blue), Af = 16 (purple). In these simulations the inflaton 
potential is quartic with A = 10 -14 and g 2 /A = 120. 
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FIG. 5: The variances of one \ a for each simulation as a func- 
tion of conformal time, r, during the stage when preheating 
becomes inefficient. The color coding is: Af — 1 (red), Af — 2 
(orange), Af = 4 (green), Af = 8 (blue), Af = 16 (purple). 



10 -is , 




lxlO 7 2xl0 7 5xl0 7 1x10 s 2x10 s 5xl0 f 
Frequency (Hz) 



FIG. 6: The present day gravitational wave spectrum from 
preheating. The color coding is: Af — 1 (red), Af — 2 (or- 
ange), Af = 4 (green), Af — 8 (blue), Af — 16 (purple), 
Af = 32 (black). In these simulations the inflaton potential is 
quartic with A = lCT 14 and g 2 /\ = 120. 



B. Quadratic inflation 

Lastly we investigate quadratic inflation. The poten- 
tial is given by 

V{4>) = \m 2 4>\ (28) 

where m — 10 _6 m p i and g 2 /m 2 — 2.5 x 10 5 . The stages 
of preheating are identical to those discussed in the pre- 
vious subsection. At the end of inflation the inflaton, 
cj>, is coherently oscillating and the modes of the matter 
fields, Xa, are subject to (approximate) Mathieu equa- 
tions. During this stage the Xa are parametrically am- 
plified by the coherent oscillations of the inflaton, and we 
expect the gravitational radiation produced to be propor- 
tional to the number of fields. This can be seen in Fig. [7] 
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FIG. 7: The present day gravitational wave spectrum pro- 
duced during the stage when the <j> field is coherently oscillat- 
ing. The color coding is: Af = 1 (red), Af = 2 (orange), Af = 4 
(green), Af = 8 (blue), Af = 16 (purple), Af = 32 (black). 
In these simulations the inflation potential is quadratic with 
m = 10~ 6 m p i and g 2 /m 2 = 2.5 x 10 5 . 

As preheating continues, the back-reaction on <f> be- 
comes important and preheating becomes inefficient. The 
A/'-dependent amplification of gravitational radiation is 
counteracted by this back-reaction and models with more 
fields cease to be more efficient. By the time the 
source term of Eq. (|15p vanishes the spectrum of grav- 
itational radiation is virtually independent of Af, see 
Fig. [H] Namely, the peak frequency and amplitude re- 
main roughly unchanged. 

C. Af and g 

The form of the potential Eq. ^ allows for the field 
redefinition, 

a 

along with a set of angles, 8 a , that cover the (Af — 1)- 
sphere. Using this set of fields, the inflation only couples 
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FIG. 8: The present day gravitational wave spectrum from 
preheating. The color coding is: TV = 1 (red), TV = 2 (or- 
ange), TV = 4 (green), TV = 8 (blue), TV — 16 (purple), 
TV = 32 (black). In these simulations the inflation potential 
is quadratic with m = 10 _6 m p i and g 2 /m 2 = 2.5 x 10 5 . 

to the norm of the field, £, 

mo = ^) + ^We 3 . (so) 

This can lead one to believe that the addition of scalar 
degrees of freedom is equivalent to changing the coupling 
in the case of a single field. The analogy is not exact, 
though, since the kinetic term must also be modified. 
Consider, for example, the case of AT — 2. In this case 
we only have one angle 9 and the kinetic term becomes, 

r - iv^w + iv^v^ + Ie 2 v^v^. (31) 

In other words, if we wish to draw an analogy with single- 
scalar dynamics, we need to absorb the angular kinetic 
term into an effective potential, 

V eS (^, = V(tt>) + \ [flW - i 2 - (32) 

Using this approach, we can see that £ has a time- 
dependent mass term, and will generically not vanish 
when vanishes. We can only see one last way to force 
our analogy to be exact; namely to set 9 = 0. Outside of 
this limiting case, 9 is a a dynamical degree of freedom 
and influences £. 

IV. CONCLUSIONS 

In the early universe, we expect the infiaton to be cou- 
pled to many of the fields present at that time. We have 
investigated the effect of these fields on the dynamics 
of preheating, and explored changes to the gravitational 



wave signature compared to the case with just one field- 
the only case investigated in the literature so far. 

Our simulations show that parametric resonance is ro- 
bust to the presence of many scalar fields coupled to the 
infiaton. In addition, we find significant gravitational 
wave production regardless of the number of fields. 

In the case of many scalars we see two subtle, TV- 
dependent effects. In both quartic and quadratic in- 
flation models we observe an enhancement in the low- 
frequency content of the gravitational wave spectrum. 
This enhancement is caused by the presence of multiple 
Xa fields, and is sourced when the <j> field is homogeneous 
and oscillating coherently. In this early phase the low fre- 
quency enhancement to the gravitational wave spectrum 
is proportional to TV. This is the stage when parametric 
resonance is most efficient. At later times the presence 
of extra matter fields induces significant back-reaction on 
the infiaton and causes parametric resonance to lose effi- 
ciency relative to the TV = 1 case. In the quartic inflation 
model this effect suppresses the high frequency content 
of the gravitational wave spectrum. The combination 
of the low frequency enhancement and high frequency 
suppression results in a slight flattening of the spectrum 
when the number of matter fields is large. Conversely, 
the effect of the many matter fields on the gravitational 
wave spectrum in quadratic inflation is significantly less 
pronounced. 

In [l3| the authors claim that the detection of a pre- 
heating signal is complementary to the detection of pri- 
mordial gravitational radiation from inflation, and that 
the energy scale of inflation can also be determined 
from the gravitational wave signature of preheating. 
The results of the present work further strengthen the 
claim that the gravitational wave signature is determined 
largely by the energy scale of inflation and only weakly 
depends on the number of matter fields that are coupled 
to the infiaton. 
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